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Abstract 

In this paper we study dynamical properties of the area preserving Henon map, as 
a discrete version of open Hamiltonian systems, that can exhibit chaotic scattering. 
Exploiting its geometric properties we locate the exit and entry sets, i.e. regions 
through which any forward, respectively backward, unbounded orbit escapes to infinity. 
In order to get the boundaries of these sets we prove that the right branch of the 
unstable manifold of the hyperbolic fixed point is the graph of a function, which is the 
uniform limit of a sequence of functions whose graphs are arcs of the symmetry lines 
of the Henon map, as a reversible map. 

1 Introduction 

The area preserving maps of the plane are discrete versions of open Hamiltonian systems 
that can exhibit chaotic scattering |^. Their phase space is an open and unbounded set, that 
is the basic assumption of the Poincare-recurrence theorem concerning their dynamics, 
is violated. Hence it is natural to wonder which points of the plane are non-wandering, and 
through which regions do some orbits escape to infinity. 



The simplest area preserving map of the plane is the quadratic Henon map 

{x,y) H-> + 1 - aa;^a;), 



(1) 



where a is a parameter. Despite its simplicity, the area preserving Henon map exhibits a very 
complex dynamics, described in a long list of papers appeared during the thirty years from 
its definition. Even last five years important new results concerning the dynamic behaviour 



of this map were reported. In |T2[, for example, Meiss presents a theoretical support for 
the study of transport through the resonance zone associated to the hyperbolic fixed point 
of the Henon map, as well as, for a numerical tool designed to compute the measure of the 
bounded orbits of the Henon map. Homoclinic bifurcations for this map are discussed in 
Tl], and the computation of its periodic orbits using anti-integrable limit in |]T3[. A method 



for computation of the self-rotation number for its orbits is given in 0]. 

In this paper we study the Henon map from the point of view of open conservative 
systems, i.e. conservative systems whose phase space is unbounded. We prove the existence 
of an unbounded forward (backward) invariant set, called exit set {entry set) through which 
any forward (backward) unbounded orbit escapes to infinity. These sets are very useful in 
the study of transport properties of Henon maps ||TT], |^, as well as of chaotic scattering in 
their dynamics |0, |10 . 



In order to locate these regions for the Henon map, we study its geometric properties, 
and prove that the right branch of the unstable manifold of the hyperbolic fixed point is the 
graph of a function, which is the uniform limit of a sequence of functions whose graphs are 
arcs of the symmetry lines of the Henon map, as a reversible map. 



2 Regions of different dynamical behaviour for Henon 
maps 

Some dynamical properties of polynomial mappings with polynomial inverse are studied in 
0. These maps are called Cremona maps and they form a group G. Planar diffeomorphisms 
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gs & G oi the form 

{x,y) ^ {y,p{y) - 6x), (2) 

where p is a polynomial function of degree at least two, and 5 7^ is the constant Jacobian 
determinant of the map gs, are called in |^ generalized Henon transformations. 

For 5 = 1 in we get polynomial area preserving maps with polynomial inverse. 

By any Cremona map F of prime degree d is affinely conjugate either to a generalized 
Henon map or to an elementary transformation e{x,y) = {ax + p{y),by + c), ab ^ 0. In 
the first case the generalized Henon map can have the polynomial p in the normal form 
p{x) = ±x'^ + terms of degree < — 2. 

For our purposes it is more appropriate to study the class of area preserving generalized 
Henon maps of the form g^^, i.e. maps defined by: 

{x,y) ^ {-y +p{x,fi),x), (3) 

where ci = 2, 3 is the degree of the polynomial family p{x, fi). Because in this case d is a prime 
number, a particular choice of the polynomial family p{x, /i) does not affect the generality 
of results. Namely, we study the quadratic Henon map H2{x,y) = {—y + fix + and 
in a forthcoming paper the cubic Henon maps H^(x,y) = {—y ± + yux + z/, x), as open 
systems. 

In our study we exploit the reversibility property of the Henon map. The area preserving 
maps (H) are reversible with respect to the involution R : ^ M^, R{x,y) = {y,x), i.e. 
H^^ = RoHdoR. Hence the map Hd factorizes as = loR, where I{x, y) = {—x+p{y, fi),y) 
is also an involution. 

Recall some properties of the reversible systems we use in our approach (a survey on the 
dynamics of reversible systems is presented in 0). 

The relation H^^ = Ro Hd o R ensures that H^^ = Ro o R also holds, for any n G Z. 
Hence is i?-reversible, too. Denote by /„ the involution defined by In = o n G Z, 
and by r„ its fixed point sets, called n-symmetry line: 

Tn = {ix,y)\iH^oR){x,y) = ix,y)} (4) 
3 



If i? is a C^-orientation reversing involution, then Fq is nonempty, and moreover it is a 
C^-curve in having no self-intersections. 

The symmetry lines are transformed by into other symmetry lines in the following 
way: 

H^iTk) = T2n+k, yn,keZ (5) 
For the Henon map (0), R and I are orientation reversing involutions and their fixed point 
sets are, respectively Tq : x = y, Ti : x = p{y,fi)/2. An /^-invariant periodic orbit of the 
map Hd is called symmetric orbit. The symmetric fixed points of the map lie at the 
intersection of the two basic symmetry lines Fq and Fi. 

Henon maps being discrete versions of open Hamiltonian systems, it is very important 
to locate the non-wandering set for such a map, and the region through which forward 
unbounded orbits escape to infinity. A point 2; G is non-wandering for the Henon map 
Hd, if for any neighbourhood V of z, there exists an integer n > 0, such that Hd{V) (IV ^ ^. 
In other words, after some time n > 0, the orbit of a non-wandering point z & V, comes 
back into the neighbourhood V. 

The set of non-wandering points of the Henon maps H^, denoted Q^Hd), is proved to be 
contained in a particular region of the phase space 0]. In this paper we locate the entry and 
the exit set for the quadratic Henon map. 

If a generalized Henon map has no fixed points, then every orbit is unbounded [Q. In 
the case when (^) has fixed points, then there is a box B, such that any orbit outside it is 
either forward unbounded or backward unbounded or both. This box is: 

B = {{x,y)\ |x|<M,||/|<M}, (6) 

where M is the largest of the absolute values of the roots of the equation \p{x, fi) \ — 2\x\ = 0. 
Hence bounded orbits lie inside the box. 

The quadratic family H2{x,y) = {—y + fix + has two i?-symmetric fixed points: 

Ze = (0, 0) and Zh = {2 — ji,2 — ji). Denote Xh '■= 2 — ji. 

We study the maps corresponding to the parameter /i G (—2, 2). In this range the origin 
is elliptic, while the second fixed point is hyperbolic. The value M in the definition of the box 
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B is, in this case M = Xh- The bounded orbits of the quadratic Henon map, corresponding 
to yU G (—2,2), are periodic orbits, all orbits filling densely invariant circles surrounding the 
elliptic fixed point whose multipliers are A, A = e^'^'"^'^, with uj ^ 1/3, 1/4, as well as orbits 
on the cantori (former invariant circles, which have broken down). 

Beside the orbits outside the box B, it is possible to exist unbounded orbits which enter it, 
spend some time in B, and then are scattered to infinity. Chaotic scattering in the dynamics 
of area preserving maps of the plane was revealed for the quadratic map f{x, y) = (x, y) [ p!0| : 

X = [x-{x + yf/4]a 
y = [y + {x + yf/A]a'^ 

where a > 1 is a parameter. 

Because any quadratic area preserving map of the plane, with nontrivial dynamics, is 
affinely conjugated to the Henon map H2-, it results from the study of the map (^, that 
in some range of the parameter /i of the Henon map can also exist a subset A in the box 
B, such that almost every point in A belongs to a trajectory that eventually leaves B (it 
is scattered to infinity). The stable set S{A) (unstable set U{A)) is the set of points in A 
whose forward (backward) orbit stays in A. The set C = S{A)r\U{A) is iJ2^invariant. If this 
set is a Cantor set which contains a dense orbit in C, of positive Lyapunov exponent, then 
C is called chaotic saddle. Hence the orbits of points on a chaotic saddle are also bounded. 
A \ C is called scattering region. A trajectory entering the scattering region typically spends 
a finite amount of time near the chaotic saddle, and then exit it, as well as the box B, 
being scattered to infinity. The scattering process is called chaotic if the motion near the 
chaotic saddle is chaotic. Chaotic scattering is a manifestation of transient chaos in open 
Hamiltonian systems. In order to quantify the scattering process in the dynamics of area 
preserving maps one defines the exit-time function, and the time-delay function The 
exit-time function associates to an initial position in A the minimal time needed to reach the 
exit set, while the time-delay function associates to a point in the entry set the minimal time 
to reach the exit set. Each of these functions gives information on the time that the moving 
particle spends near the chaotic saddle, bouncing between its points. Because the points of 
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the chaotic saddle do not leave A, the exit-time function has a Cantor set of singularities. 

As far as we know, the entry and exit sets of the polynomial area preserving maps, 
analyzed from the point of view of chaotic scattering, have not been located, and as a 
consequence the exit-time of a scattered particle zq was associated in a subjective way, 
namely as the the minimal time n such that the n*^ iterate of the map applied to Zq is 
outside a large disk centered at the origin |jlO|, ^. Moreover the time-delay function could 



not be evaluated because the entry set was unknown. 

Next we show that for the quadratic Henon map we can determine analytically an exit 
set (entry set) of the system, i.e. an unbounded forward (backward) invariant set E [F), 
such that any forward (backward) unbounded orbit enter E {F) , and as a consequence goes 
to infinity through this set. 

The entry and exit sets are also useful in the study of transport through the resonance 
region, bounded by segments of stable and unstable manifolds of the hyperbolic fixed point 
of the Henon map [0, ^ . 



3 The exit and entry set of the quadratic Henon map 

In order to detect the exit and entry sets of the Henon map we study its geometric properties. 
For simplicity we drop the index 2 from its name H2. 

Proposition 3.1 The Henon map H{x,y) = {—y + /ix + x^,x), /i G (—2,2), has an un- 
bounded forward invariant set Eq C {(x, y) G M^l x > 0,?/ > 0}, i.e. H{Eq) C Eq. Its 
boundary OEq is CU Fq, where C is the curve parameterized by (X(a), Y{a)), a G (0, 1], with 

- - a/i + 1 

A a = 

^ ^ a (8) 

Y{a) = aX{a), 

and F[f = {{x,y) \ y = x, x > Xh}. 
Proof 
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Let L{a) be the infinite radius (semi-line) y = ax, starting from tlie elliptic fixed point 

(0,0). Its image under H is an arc of parabola, x = y'^ + {fi — a)y. Only the radii 

Li(a) defined by x > 0, a > 0, intersect their images H{Li{a)) at the points /(a) = 
a"^ — afi + 1 — afi + V 



a 



(Fig-i- 



The point I{a) is the if-image of the point P{a) = I , a — a/i + 1 . For any 



a — a/i + 1 



a 



H G (—2,2), /(a), with a > 0, belongs to the set S = {{x,y) | x > 0, y > 0}, as well as the 
sub-arc of the parabola x = + (/i — a)y corresponding to y > yi(a)- This suggests to look 
for a forward invariant set in S. 

Any point Q{xo, axo) G Li{a), < a < 1, with xq > Xp(^a) is mapped by H to the point 
(a;i,a'a;i) G Li{a') (see Fig.|l]), with 



a'{xo, a) = -^^ < a (9) 

Xo + yu - a 



Let us show that the set points 

Eo = {ix,ax)\x> ^ ~"^^^ ,0<a<l} (10) 

is forward invariant. 

a'^ — a' a + 11 

For, take a pomt [Xo,axoj m Eq. By [m, = \- Xq — a. From 

a' xo + /i — a 

a'^ — a'fi +1 — a/i + 1 

a < a, we have ; < Xq. For < a < 1, > 1 — /i + a. Hence 

a a 

Xq > 1 — u + a, I.e. xo + — a > 1, and thus xi > xq or equivalently xi > . 

a' 

This means that (xi,a'xi) = H{xo,axo) G -Eq, for any (xo,axo) G -Eq- □ 

Remark 3.1 In canonical polar coordinates, a point (ro, ^o) ^ Eq is mapped by H to a point 
{ri,6i), with 6*1 < Oq. Any point (ro, ^^o) ^ (0, cxo) x [7r/2,27r], zs mapped to a point (ri,6'i), 
where 9i is obtained from 9q by a positive translation (rotation) on the circle. 

Proposition 3.2 The forward orbit {xn,yn) = H"'{xo,yo), n G N, of a point (xo,?/o) ^ Eq 
has the property that both sequences (x„), (y„) are increasing and unbounded. 

Proof 
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Projection onto the first factor of the set Eq is {xh, oo). Hence Xq > Xh- From the proof 



of the Proposition ^]T] it results that the sequence (a;„) is increasing. Moreover 

Xn+l+Xn-l=p{Xn), (11) 

with p{x) = fix + x"^. The equation (|l^) is the Henon map expressed as the second order 
difference equation. The sequence (x„) is unbounded, because otherwise it converges to a 
point X* satisfying, by (pA|), the fixed point equation associated to the Henon map, p{x*) = 
2x*, i.e. X* ^ pT^{Eo). 

In a similar way one proves that the sequence is increasing and unbounded 
satisfies the same second order difference equation). □ 

As a consequence, in canonical polar coordinates, the orbit (r„,6'„), n G N, of a point 
(ro, 6*0) G Eq, has r„ = (x^ + y^)/2 oo, i.e. the orbit goes to infinity. 

Let A" > 1 be the expanding multiplier of the hyperbolic fixed point z^, and v = (A", 1) a 
corresponding eigenvector. The branch of the unstable manifold W^{zh) having the direction 
and the sense of v for the tangent semi- line / at Zh is denoted by W^{zh), and called the 
right branch of the unstable manifold. 

Proposition 3.3 The right branch W!^{zh) of the unstable manifold of the hyperbolic fixed 
point Zh is included in the set Eq. 

Proof 

A simple computation shows that the segment I = {z E M.'^\z = z^ + tv , < t < 1} of the 
tangent at Zh, to the unstable manifold, is included in the set Eq. Hence, there is a 5 > such 
that an arc of the local unstable manifold Wl^^{zh) = {z E M?\ H~'^{z) G Bs{zh),^n > 0} 
is included in Eq {Bs{zh) is the (5-ball centered at Zh). Denote by L" this arc. The right 
branch of unstable manifold is then W^{zh) = |J H"'{L^), and by the forward invariance of 

n>0 

Eq, it is included in this set. □ 



The Henon map H being reversible, its inverse is conjugated with H by the involution 
R{x,y) = {y,x). Hence the set Fq = R{Eq) is a backward invariant subset, i.e. H^^{Fq) C 
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Fq. Its boundary is OFq = C U Fq, where C = R{C). Furthermore, the right branch, 
W^{zh) = R{W^{zh)), of the stable manifold W%Zh), is included in the set Fq. 

Next we determine the maximal forward, respectively backward, invariant set in S* = 
{{x,y) \ x > 0,y > 0}. For, we exploit the fact that Tq = Eq n Fq, i.e. orbits starting 
from points Zq = (a;o,a;o), with xq > Xh, are unbounded in both directions. Moreover, 
F'i^n := H-^{r^Q) C Fo, and F^„ := H^iT^) C Eq, V n G N. 

Because Fq C Eq it follows that the sequence of sets (En), n E N, with dEn = CU F'^ig^, 
is a sequence of forward if-invariant sets. We prove that it is an increasing sequence and 
its limit has as a partial boundary the right branch, W^{zh), of the stable manifold of the 
hyperbolic fixed point Zh- 

Proposition 3.4 The right branch, W^{zh), of the stable manifold is the graph of an in- 
creasing C^-function ip : [x/^, oo) [xh,oo), which is the uniform limit of a sequence of 
increasing -functions ipn '■ [xh, oo) [xh, oo), n G N, and: 
i) Vn{xh) = Xh,y ne N; 
a) ipn{x) > iprn{x), W u > m, and X > Xh; 

Hi) the graph of the function ipn is the partial boundary T^2n forward invariant set En- 

Proof 

Adapting a result from p[ to the case of area preserving Henon map, we have that the 
right branch W^{zh) of the unstable manifold is the graph of an increasing C^-function 
g : [xh,oo) {xh,oo). The right branch of stable manifold W^{zh) = R(W^{zh) being the 
symmetric of W^{zh) with respect to the line y = x, is the graph of the function (f = g^^, 
and this function is also increasing. 

The partial boundary Fq of the set £"0 is the graph of the function (po = (/^q ^, (fo{x) = x, 
X > Xh- Its image under H is the set of points H{x,ipQ^{x)) = {—ipQ-^^x) +p{x),x), with 
X > Xh, and p{x) = fxx + x^- Denote by ^i{x) := p{x) — v^q ^(x). Obviously, ^i{xh) = Xh, 
and ipo{x) < ipi{x), V x > Xh- The derivative of ipi has the property: 

if[{x) > 1, \/ X > Xh (12) 
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Hence it is an increasing, that is invertible function, and the graph of its inverse is while 
the graph of is T^^j = 

is thus the if-image of the graph of the function yj]"^, i.e. the set of points H{x, ip^-^^x)) 
= {p{x) — Lpi^{x),x), X > Xh- Denote by the function defined by '^2{x) := p{x) — Lp^^{x). 
^2ixh) = Xh- Because ^po{x) < ^i{x), V x > Xh, is equivalent to ^i^{x) < lPq^{x), V x > Xh, 
we get from the definition of the functions ipi, i = 1,2, that f2{x) — fi{x) = ipQ^ — ipi^ > 0. 
The property (|1^) and the fact p'{x) > 2, for any x > Xh, ensures that f'2{x) > 1, Vx > x/^. 
Hence the graph of ip2'^ is T'l, and the graph of is T'^^. By induction, we get that the 
partial boundary T'l2n of the set En, is the graph of an increasing C^-function : [xh, oo) 
[xh, oo), defined by: 

ipn{x) := p{x) - (fn-i, V X > Xft, (13) 

and it fulfills the conditions: (pn{xh) = Xh, ^Pn{x) > ipn-i{x), V x > Xh, and v^'„(x) > 1. 

No arc r^2fc '^^^ intersect the stable manifold W^{zh) in Fq. To prove this property, 
suppose that there exists a point z G W^{zh) fl r'J:2fc, z ^ z^. Hence lim„^oo -^^"(-z) = Zh- 
Because Tt2k ^ r_2fe, by (||) we have [H^'^^ o R){z) = z. At the same time, Rz e Wizh), 
and lim„^oo -^^^"(-R^) = or equivalently lim„^oo = Zh- Hence the 
orbit of 2; is a homoclinic orbit to Zh ^ Fq, which contradicts the fact that 2; is a point in the 
backward invariant set Fq. 

Because W^{zh) cannot intersect any arc of symmetry line r_2n in -Fq, we have that 
(Pn{x) < (p{x), for any x > Xh- Hence the sequence of functions converges uniformly 
to an increasing C^-function : [xh, 00) —>■ [x^, 00). By (|T^ this function has the property 
that p{x) — ip{x) = ip'^^x). Hence a point {xo,4'{xo)) of its graph is mapped by the Henon 
map to H{xo,4'{xo)) = (p(xo) — ■ip{xo),Xo) = {ip~^{xo),Xo). Denoting xi = ip~^{xo), we have 
H{xo,4'{xo)) = {xi,4'{xi)), i.e. the graph of the function ip is forward invariant under the 
action of the map H. If (x„, ip{xn)) = H^{xo, i/j^xq)), then (x„_i, '?/'(x„_i)) = H~^{xn, ipi^Xn))- 
By Proposition p.2| and the reversibility of the system we get that x„_i > x„, V n > 0, and 
Xn Xh- Hence if"(xo, ^'(^^o)) {xh,Xh), as — * 00, that is, the graph of ip is the right 
branch of the stable manifold of the fixed point Zh, i.e. ip = if - a 
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Corollary 3.1 The sequence of forward invariant sets En having the boundary dEn = CU 
r^2n fl?^ increasing sequence and its limit E = Un>oEn is forward invariant, too, and has 
the boundary dE = £ U W^{zs). 

The set F = R{E) is backward invariant, and OF = U W^{zs)- 

The set of points G = E \^ F = {{x,y) \ x > Xh,^p~^{x) < y < '^{x)} is H -invariant, i.e. 
H{G) = G. 

Because the orbit of a point in G is unbounded in both directions, the exit and entry sets of 
dynamical interest for transport processes, as well as for chaotic scattering are E' = E\G, 
respectively F' = F \ G. The boundary for E' is dE' = C U W^{zh), while for F' is 
OF' = C'UW^izh). 

The position of these sets in the phase space of the map is illustrated in Fig.^, and 
the position of the invariant manifolds of the hyperbolic fixed point, with respect to these 
regions is shown in Fig.^. 

By straightforward computation we get that the preimage of the curve C, as well as 
the image of intersects the box B. Hence there exist in the box B, points z, such 
that t^{z) = min{n > | H^{z) G E'} = 1, respectively points z' with t~{z') = min{n > 
I H^"'{z) G F'} = 1. This means that indeed E', F' are exit, respectively entry, set for the 
Henon map, as a map having a scattering region. 

Knowing these sets, a typical scattering experiment consists in launching particles z E F' 
and observing their trajectory. If such a trajectory enters the box B and then exits and 
reaches i?', the time-delay is the time it takes to a particle to go from the launching position 
until it reaches the exit set E. To a particle z E F, which do not enter B, one associates 
the time 0. The existence of particles in F' that do not enter B is ensured by the fact that 
H{C') intersect the complementary set of B, not only the set B. 

In Fig.|^ is illustrated the time-delay function of points in the entry set F' belonging to 
the segment y = 3.61, x G [2.66,2.74]. 

Note that the points of the curve £ are mapped to points in E'. Namely, the points on C 
are the points P{a) in the Proposition and they are applied to points /(a) G E', i.e. to 
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points on the same semiline of slope a (see Fig.|T]). Moreover, a simple analytical computation 
shows that points on the semi-line ?/ = 0, x > Xh-, as well as points in the complement of 
the box B, lying between this semi-line and the curve £, are applied in few iterations of the 
map into the set E' . Taking also into account the Remark |3]^, we conclude that any point 
from the exterior of the set having a forward unbounded orbit, reaches the exit set in 
a small number of iterations. 

4 Conclusions 

Studying geometric properties of the quadratic area preserving Henon map, we have located 
in the phase space its exit and entry sets. The knowledge of these sets allows a more rigorous 
study of transport properties, as well as of chaotic scattering processes in the dynamics of this 
map. It appears that unlike classical Hamiltonian systems, exhibiting chaotic scattering, in 
the dynamics of the Henon map, the particles cannot enter the scatterer from any direction 
of the phase space, and cannot exit from the system in any direction. The entry (exit) set is 
located in a neighbourhood of a branch of stable (unstable) manifold of the hyperbolic fixed 
point of the map. 

In a forthcoming paper we determine these sets for generalized Henon maps of degree 
three, which have either an elliptic fixed point and two hyperbolic fixed points, or a single 
elliptic point, but a two periodic hyperbolic orbit. 
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Fig. 1: The action of the map if on a semi-hne Li{a), a <1. 
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Fig. 2: a) Regions of different dynamical behaviour of the Henon map, corresponding to 
ji — —0.78; b) The invariant manifolds of the hyperbolic fixed point. Within the sets E' [F') 
one can see foldings of the left branch of the unstable (stable) manifold of the point z^,. 
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Fig. 3: The time-delay function defined on the segment y — 3.61, x e [2.66,2.74] included 
in the entry set F' of the Henon map corresponding to /i — —0.75. 
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